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Abstract: We present a study of loop-level lepton flavor vioalting signals in models with 
' heavy Majorana neutrinos and in the supersymmetric extension of the standard model. 

The attention is focused to the e~e~ option of the next generation of linear colliders and 
■ its potential of discovering new physics is emphasized. 



particle interactions. Besides confirming and allowing precision measurements on new 
^ . physics which, hopefully, will be discovered at the LHC, they offer by themselves new 

physics opportunities thanks to the possibility of 77, 7e~, e~e - collisions. These options 
will allow to test, with higher sensitivity couplings of the standard model (and alternative 
theories) which cannot be studied in hadronic or e + e~ collisions. Moreover the possibility 
of employing electron beams with a high degree of longitudinal polarization will be essential 
to enhance signals of new physics. 

The e~e~ mode with lepton number L = +2 of the initial state is particularly suitable 
for studying lepton and flavor number violating reactions (LFV). The "merits" of this 
option have been recently recalled by the ECFA/DESY working group in Ref. JlJ. Among 
these there is (1) the possibility to identify Majorana neutrinos of masses in the TeV range 
through the "inverse neutrinoless double beta decay" reaction e~e~ — > W~W~, and (2) 
the sharper onset of slepton production threshold respect to the e + e~ option. Here we 
present a study of two loop-level reactions strictly connected with these two "merits" : we 
first consider seesaw type models with heavy Majorana neutrinos at the TeV scale and 
study the reactions e~e~ — ► £~l~, and then similar reactions e~e~ — ► £~e~ (£ = fi,r) 

'Speaker. 





A ■ 1. Introduction 

o 

i-G ■ Linear e + e colliders are very important step towards the understanding of high energy 
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in supersymmetric models where LFV is due to slepton mixing. The importance of the 
e - e~ mode for studies of LFV due to slepton mixing and the advantages of the threshold 
behavior was first noted by Feng in various papers ||, ^ |], [B|, [|. 

We briefly discuss the standard model background. More details on the calculations 
and numerical tools used can be found in Refs. |7|. |8[]. 

2. e~e~ — > t~t~ (£ = fi, r) through heavy Majorana neutrinos 

To give a detectable signal, heavy Majorana neutrinos (HMN), besides having masses in the 
TeV range, must have interactions which are not suppressed by the mixing matrices with 
light states as instead happens in the one family seesaw mechanism, where 9 ~ ^Jm v jM^ . 
With three generations, more free parameters are at our disposal, and the "two miracles" 
of not so large masses and non negligible mixing, are obtained imposing suitable relations 
among the elements of the matrices mo and Mr: examples of these models were proposed 
some time ago in Refs. |, |fi], [lit] and in the 

more recent paper, Ref. [{T2; 1, whose authors 
suggest "neutrissimos" as the correct name for these particles. According to other authors, 
Ref. these models are based on fine-tuned relations, but are shown to be not in 

contradiction with any experimental bound. We do not enter in such a theoretical dispute 
and assume "neutrissimos" in the TeV range and study the phenomenological consequences. 

Experimentally one cannot put bounds on the single mixing matrix elements, but 
on some combinations of them, assuming that each charged lepton couples only to one 
heavy neutrino with significant strength. Light-heavy mixing has to be inferred from low- 
energy phenomenology and from global fits performed on LEP data identifying the following 
effective mixing angles s\ = • | -S^iTVj 1 2 = sin 2 6*^ £ with upper bounds 

s 2 e < 0.0054, si < 0.005, s 2 T < 0.016, (2.1) 

where B is the mixing matrix appearing in the charged current weak interaction lagrangian. 
Under these assumptions, the coupling of neutrissimos to gauge bosons and leptons is 
numerically fixed to gB^ v where g is the SU(2) gauge coupling of the SM. Since the 
width of the heavy states grows as Mjy, at a certain value it will happen that r^r > Mjsj, 
signaling a breakdown of perturbation theory. The perturbative limit on Mn is thereby 
estimated requiring < Mjv/2, which gives an upper bound of ~ 3 TeV |1| for the 
numerical values given in Eq. fl2.ip . 

We study the process e~e~ — » (£ = fi,r), which violates the L e and lepton 

numbers. The amplitudes are computed in the 't Hooft-Feynman gauge where there are 
graphs with WW, <p<p and <pW exchange, 4> being the Goldstone boson, as shown in Fig. |l|. 
The numerical computation of the four-point functions was performed using the loop- 
tools [ 16 1 package. A similar study was done in Ref. [17] using the approximation where 
all external momenta in the loops are neglected relative to the heavy masses of the gauge 
bosons and Majorana neutrinos, that is qi = k + k{ ~ k, where k is the integration variable 
in the loop integrals and ki are the sums of external momenta appearing in the propaga- 
tors. This approximation for the four-point functions is good at low energies, such as in 
decay processes of heavy mesons, or when << M, M being the highest mass running 
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Figure 1: In (a-d) the Feynman diagrams, 
in the 't Hooft-Feynman gauge, contributing 
to e~ eT — > l~t~ (I = /j,t) via heavy Ma- 
jorana neutrinos. <fi is the unphysical Gold- 
stone boson. 



Figure 2: (a) the ratio R a is plotted as a 
function of y/s. Solid line, Mn 4 — 
lliii GeV; dotted, M Ni = M Nj = 3 TeV. 
(h) Solid line, our calculation, dashed line 
that of Rcf. 0. 



in the loop. In this way the cross section presents a linear growth with s which breaks 
unitarity: therefore, in order to make quantitative predictions with the correct high energy 
behavior, full dependence on the external momenta of the four-point functions has to be 
considered. Theoretically, according to the 'Cutkosky rule', one expects an enhancement at 
y/s — 161 GeV ~ 2M\y, the threshold for on-shell WW gauge boson production, at which 
the four-point functions develop an imaginary part. In Fig. Q(a) the ratio R a = crtot/o~o of 
the integrated total cross section a tot to uq, the cross section of the low energy calculation 



of Ref. [17 1, is plotted for sample values of the Major ana masses. The enhancement due 
to the threshold singularity of the loop amplitude is more pronounced for values of Ma- 
jorana masses close to M\y and is drastically reduced increasing Mjv 4 ~ Mjy^ to O(TeV). 
As R a — ► 1 at \/s « M\\r in all the cases, the agreement of our full calculation with the 



result of Ref. [17] in the regime of low energies is evident. The threshold effect appears to 
be quite spectacular only for values of Majorana masses which correspond to cross sections 
too small to be measured even at a next linear collider. In Fig. §(b) the effect of the 
threshold singularity in the loop integral is shown reporting absolute cross sections for a 
particular choice of Majorana masses: = 150 GeV and Mjv- = 450 GeV. The low 

energy approximation (dashed line), obtained neglecting external momenta in the loop, is 
inadequate when the energy of the reaction increases to values comparable with the masses. 
Increasing the energy, after reaching a maximum, the cross section starts to decrease until 
the asymptotic behaviour 0(1/ s 2 ) of the loop integral is reached. This happens for every 
value of heavy Majorana neutrino masses and we checked numerically that, as expected, for 



M 

> 
EC 

o 
o 



- 3 - 



International Workshop on Astroparticles and High Energy Physics Mirco Cannoni 



o(fb) 
10° 



(a) 



(b) 



i 1 1 inn — 1 1 iiiiij — it™ - i 1 1 iiiiij — 1 1 iiiiij — i mm 



M r ITeV,M j =3TeV wl 
M.=1TeV,IMTeVn 




/ __ M=M=3TeV(ln ; 



I I 




I l/lllllll I 



0.2 



0.15 



CD 

S 0.1 
o 





0.05 



1 2 1 4 1 2 3 4 

10 10 10 10 10 10 10 10 

Vs(GeV) 

Figure 3: Total cross sections as function of 
y/s. (a) the solid curve referes to the case of 
e~e~ -> t~t~ with M Ni = M Nj = 3 TeV, 
while the dashed line referes to (e~e~ — > uu) 
with the same values of Majorana masses, 
(b) the Majorana masses are changed to 
somewhat lower values: Mjf t = 1 TeV, 
M Nl = 3 TeV. 
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Figure 4: Angular distribution in the polar 
angle of the outgoing lepton for different val- 
ues of the center of mass energy, y/s in the 
case of e~e~ — > t~t~ with Mjv 4 = Mm = 3 
TeV. The curves are not exactly constant, 
and using an appropriate scale they show 
small deviations from a stright line, remain- 
ing left-right symmetric. 



higher masses the asymptotic regime is reached at higher values of y/s. In fact from Fig. [| 
we note that the cross section grows with increasing HMN masses. The main contribution 
comes from the graph with two Goldstone bosons since their coupling is proportional to 
M]\f i . Moreover the chiral structure of the coupling selects the mass term in the numerator 
of the Majorana neutrino propagators. When these masses are much larger then the other 
quantities, the amplitude scales like Mj^Mjy /Mj^Mjy ~ M^.M^, i.e. is proportional 
to the square of the heavy masses. This fact is the well known non decoupling of heavy 
fermions in theories with spontaneous symmetry breaking (similarly in the SM the top 
quark gives sizable radiative corrections owing to its large mass and a quadratic non de- 
coupling). In Fig. |3](a) the cross section is plotted for masses up to the perturbative limit, 
using the maximally allowed value of the mixing. We see that for = Mm = 3 TeV the 
signal does reach the level of 10 _1 , 10 -2 fb respectively for the (rr) and the (nfi) signals at 
y/s = 500 GeV, which for an annual integrated luminosity of 100 fb -1 would correspond 
respectively to 10 and 1 event/year. At higher energies, O (TeV), one could get even larger 
event rates (30 and 3) respectively. The solid curve refers to e - e~ — > t~t~: this is largest 
because the upper limits on the mixing are less stringent. One can also see the onset of the 
asymptotic regime at -y/s 3 TeV. Fig. ||(b) shows that the cross section quickly decreases 
as lower Majorana masses are considered. 
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As even in the more optimistic cases event rates are quite modest it is important to 
check how the signal cross-section is affected by kinematic cuts on the angle of the outgoing 
leptons. The angular distributions turn out to be practically constant as shown in Fig. |||. 
They are forward-backward symmetric because both the t and u channel are present. The 
absence of a strong dependence on the polar angle is due to the fact that within the range 
of the parameters used here the contributing four point functions depend very mildly on 
the kinematic variables (u and t). This behaviour can be most easily understood using 
helicity amplitudes. In the limit of massless external particles the process is dominated 
by a well defined helicity amplitude: eiei — ► £l^l- In the center of mass frame this is a 
S-wave scattering with J z = 0, meaning that the scattered particles are emitted back to 
back but without a preferred direction relative to the collision axis (z). Thus this signal is 
characterized by practically flat angular distributions and as a result the total cross section 
is quite insensitive to angular cuts. We have verified that with | cos 6*| <0.99 the change 
in o~t is ~ 1% for all energies considered, while using | cos#|<0.95 the total cross-section 
decreases by « 5%. The reduction of the total cross section is measured almost precisely by 
the reduction of the phase space, meaning that the angular distribution is constant up to 
~ 0.1%. Thus it can be concluded that the number of events will not be drastically affected 
for any reasonable choice of experimental cuts. The possibility of employing electron beams 
with a high degree of longitudinal polarization as planned for the new colliders will be 
essential to discover rare signals: in our case only an helicity amplitude contribute, so that 
using left-polarized electron beams will single out the essential spin configuration and we 
can gain a factor ~ 4 on the previous cross sections, because the avarage on the initial 
spins is not needed. This consideration applies also in the next Section. 

3. e~e~ — > £~e~ (£ = fi, r) in R-conserving supersymmetric models 

Slepton and squark mass matrices in soft-breaking potential, for example m ~ , are in general 
complex non-diagonal matrices. The diagonalization implies the presence of generational 
mixing matrices at the lepton-slepton-gaugino vertices. These couplings can originate too 
high rates for rare unobserved flavor changing processes. The SUGRA boundary condi- 
tions and the renormalization group equations (RGE) evolution originate flavor-diagonal 
sfermion mass matrices, so besides offering a mechanism for reducing to a manageable 
number the soft parameters, also cure the 'flavor problem'. On the other hand, if one 
wants to study the phenomenology of the MSSM without referring to a particular high 
energy scenario, one can consider general mass matrices and accept their flavor violating 
entries as large as the experimental bounds allow them. 

However when the seesaw mechanism is embedded in the MSSM with mSUGRA bound- 
ary conditions at high energy and heavy Majorana neutrinos are out of the colliders reach, a 
new source of LFV arises. The seesaw mechanism requires that the superpotential contains 
three SU(2)l singlet neutrino superfields iVj with the following couplings |[^, |i~9| ]: 



W = (Y v ) kl e l3 HlN k L\ + i(M fl )y^J^. 



(3.1) 
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Here if 2 is a Higgs doublet superfield, Li are the SU(2)l doublet lepton superfields, Y v is 
a Yukawa coupling matrix and Mr is the SU(2)l singlet neutrino mass matrix. With the 
additional Yukawa couplings in Eq. ( |3.lD and a new mass scale (Mr) the RGE evolution 
of the parameters is modified. Assuming that Mr is the mass scale of heavy right-handed 
neutrinos, the RGE evolution from the GUT scale down to Mr induce off-diagonal matrix 
elements in (m~)y. At the GUT scale we assume the universal conditions m| = m~ = 
m 2 H2 = tuq and A v = amoY u , thus the correction of the off-diagonal elements are (in the 
leading-log approximation) : 

«h * ~ (3 + al)ml(Y}Y v ) ij In (^g^) , (3.2) 

which depends crucially on the non-diagonal elements of (YjY u )ij, the square of the neu- 
trino Yukawa couplings. The main point is that these elements can be large numbers 
because in the seesaw mechanism they do not directly determine the mass of the light neu- 
trino, but only through the seesaw relation m v ~ w? d /Mr = v 2 Y 2 /Mr. On the other hand 
the same effect on the mass matrix of SU(2)l singlet charged sleptons (jn^ij is smaller: 
in fact in the same leading- log approximation of Eq. (|3.2|), the corresponding RGE do not 



contain terms proportional to yJy u , since the right-handed leptons fields have only the 
Yukawa coupling Yg, which completely determines the Dirac mass of the charged leptons 
and these are known to be small numbers. Thus the off-diagonal elements (rrv^jij can 
be taken to be ~ 0. The slepton mass eigenstates are obtained diagonalizing the slepton 
mass matrices. The corresponding mixing matrices induce LFV couplings in the lepton- 
slepton-gaugino vertices $ L .ULij&LjX" The magnitude of LFV effects will depend on the 
RGE induced non diagonal entries and ultimately on the neutrino Yukawa couplings (Yv)ij- 
The rate of LFV transitions like l{ — ► ij, i 7^ j, t = e, /j>,t induced by the lepton- 
slepton-gaugino vertex is determined by the mixing matrix U^j which, as stated above, 
is model dependent. In a model independent way, however, one can take the lepton, 
slepton, gaugino vertex flavor conserving with the slepton in gauge eigenstates, so that 
LFV is given by mass insertion of non diagonal slepton propagators. In a similar spirit 
our phenomenological approach will be quite model independent and, in order to keep the 
discussion simple enough, the mixing of only two generations will be considered. Thus the 
slepton and sneutrino mass matrix are: 

(ffi ^ /\ Tti^ 1 \ 
Am 2 rh 2 ) ' (3 ' 3) 

with eigenvalues: ih\ = rh 2 ± Am 2 and maximal mixing matrix. Under these assumptions 
the LFV propagator in momentum space for a scalar line is 

~ ~a , i ( 1 1 \ Am 2 

J I \p A — p z — m_ J (p A — m+)(p z — m_) 

and the essential parameter which controls the LFV signal is 
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Figure 5: Feynman diagrams for e~e~ col- 
lisions. The full black dot in a scalar line de- 
notes again the lepton flavour violating prop- 
agator. Exchange diagrams are not shown. 



Figure 6: Differential cross sections as a 
function of the scattering angle. Values of 
the parameters: Mi = 80, Mi = 160, = 
m c = 100 GeV and Am 2 = 6000 GeV 2 . 



In this scenario we consider the reaction 

e~(pi,Ai)e~(p2,A 2 ) -> £~{p3, A 3 )e~(j>4, A 4 ). (3.6) 

whose Feynman diagrams are shown in Fig. ||. Here Xi denotes the helicity of particle 
i. The total unpolarized cross-section (averaged over initial spins) is a = (1/4) ^ - o~y 
The signal is suppressed if neutralinos and charginos x 0,ziz are Higgsino-like, since their 
coupling is proportional to the lepton masses. For the same reason left-right mixing in the 
slepton matrix is neglected. Therefore it is assumed that the two lightest neutralinos are 
pure Bino and pure Wino with masses Mi and M<i respectively, while charginos are pure 
charged Winos with mass M2, Mi and M2 being the gaugino masses in the soft breaking 
potential. Numerical results are obtained using the mSUGRA relation Mi ~ O.5M2 for 
gaugino masses while Am 2 and the slepton masses are taken to be free phenomenological 
parameters. The parameter space is scanned in order to identify the regions which may 
deliver an interesting signal. The contributing amplitudes are 

Mei = M{ele' L ^i L e L ), 
M E 2 = M{ele- R -> r L e R ), 

M E3 = M{e- R e~ L - C L e R ). (3.7) 

The corresponding differential cross sections are plotted in Fig. || Mei has J z = 0, is 
flat and forward-backward symmetric because of the antisymmetrization. M E2 and M E3 
describe P-wave scattering with J z = +1 and J z = — 1 respectively: in order to conserve 
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Figure 7: erg 1 as a function of y/s with slep- 
ton and gaugino masses as in Fig. 6. Solid 
lines: Am 2 increasing from 100 GeV 2 to 900 
GeV 2 in steps of 100. Dashed lines: from 
1000 to 8000 GeV 2 in steps of 1000. 



Figure 8: Total cross section for the ampli- 
tude El in function of Sll for y/s — 2%. 
The values of the other parameters are given 
in the legends. 



angular momentum Me2 must be peaked in the forward direction while A4e3 favours 
backward scattering. Both M.E2 and M.E3 are orders of magnitude smaller than M.ei- 
The signal cross section is to a very good approximation given by the amplitude Mei- 
Since it is almost flat the angular integration will give a factor almost exactly equal to two. 
This again shows the importance of the option of having polarized beams. If both colliding 
electrons are left-handed one singles out the dominant helicity amplitude and a factor four 
is gained in the cross section relative to the unpolarized case. This may be important in 
view of the relatively small signal cross section one is dealing with. The analysis of the total 
cross section as a function of y/s is the following (see Fig. |^): the box diagrams dominate 
at y/s = 2m l where a changes of orders of magnitude giving a sharp peak that is smeared 
only by large values of Am 2 , while penguin diagrams give a substantial contribution only 
at higher energies. This can be easily understood considering the threshold behavior of the 
cross section for slepton pair production [||, |2(J : defining = — 4m? / s the selectron 

velocity, the amplitude of the intermediate state e^e^ — * e L e~2 behaves like 0, while for 
the other two cases it goes like (3 3 . The dependence of oei on 5ll is shown in Fig. || 
With SUSY masses not much larger than ~ 200 GeV the signal is of order 0(1O -2 ) fb for 
5 LL > O(10 _1 ). In addition the cross section is practically angle independent and thus 
insensitive to angular (or tranverse momentum) cuts. 

The phenomenological points of the SUSY parameter space corresponding to gaugino 
masses (M x , M 2 ) = (80, 160) GeV or (100, 200) GeV and to slepton masses m L = 
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Figure 9: Scatter plot in the plane (SlLjItil) of: (a) the experimental bounds from /i — > ej and 
t — * /i7 (allowed regions with circular dots); (b) regions where the signal can give at least one event 
with two different values of integrated luminosity (squared dots), for two sets of gaugino masses. 
Each signal point is calculated at y/s = 2rfiL. 



100 - 200 GeV and 5 L l > 10" 1 (which implies Am 2 > 10 3 GeV 2 ) can give in the e e 
mode a detectable LFV signal (e~e~ — ► £~e~) although at the level of 0(1 — 25) events/yr 
with Lq = 100 fb _1 . Higher sensitivity to the SUSY parameter space could be obtained with 
larger Lq. It is interesting to note that this light sparticle spectrum, which is promising for 
collider discovery, is also preferred by the electroweak data fit. In Ref. [22] it is shown that 



light sneutrinos, charged left sleptons and light gauginos improve the agreement among 
the electroweak precision measurements and the lower bounds on the Higgs mass. On the 
other hand the experimental bounds on rare lepton decays jx, r —* ej set constraints on the 
LFV violating paramters Am 2 or 5ll- the upper bounds on the branching ratios define 
an allowed (and an excluded) region in the plane {SlLi m h) which are computed using the 
formulas given in Ref. [Gl| (adapted to our model) for the LFV radiative lepton decays. 
These regions have to be compared with those satisfying the "discovery" condition 

L a(5 LL ,m L )>l. (3.8) 

Such a comparison is shown in Fig. |9| from which it emerges that: (i) For the e~e~ — ► £~e~ 
process there is an observable signal in the upper left corner of the (5iL, m L) plane. The 
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extension of this region depends on Lq. (ii) The bound from r — > ej does not constrain the 
region of the (Sll^l) plane compatible with an observable LFV signal and therefore the 
reaction e~e~ — > r~e~ could produce a detectable signal whithin the highlighted regions 
of the parameter space (upper-left regions in the (Sll^l) plane). (Hi) As regards the 
constraints from the fx — * e^f decay the allowed region in the (SlLi roi) plane is shown by 
the circular dark dots (red with colour): the process e~e~ — ► fi~e~ is observable only in a 
small section of the parameter space since the allowed region from the fi — > e~f decay almost 
does not overlap with the collider "discovery" region except for a very small fraction in the 
case of gaugino masses (Mi = 80 GeV and M2 = 160 GeV). The compatibility of values 
of 8ll ~ 1 is due to a cancellation among the diagrams that describe the i — ► £'7 decay in 
particular points of the parameter space. 



4. Standard model background 

These signals have the unique characteristic of a back to back high energy lepton pair and 



no missing energy. Sources of background were qualitatively discussed in Ref. [11]. Here we 
discuss the reaction e~e~ — > h , e u e W~*W~* followed by the decays W~W~ — > l~vit~ % , 
with four neutrinos and a like sign-dilepton pair that can be of the same or different flavour. 
This appears to be the most dangerous background, as it produces two leptons and missing 
energy, and therefore it is analyzed in more detail. 

Figure pi] shows the total cross section e~e~ — > v e v e W~W~ calculated with the Com- 
pHEP package plj| , that allows to compute numerically the 17 Feynman diagrams con- 
tributing at tree level. Above the threshold for W~W~ gauge boson production the cross 
section rises rapidly by orders of magnitude, becoming almost constant at high energies. 
In the region y/s ~ 250 — 400 GeV it increases from 10~ 2 fb to 1 fb. In order to get 
an estimate of the cross section for the six particle final state process, the cross section 
o~(e~e~ — * W~W~vv) has to be multiplied by the branching ratio of the leptonic decays of 
the two gauge bosons, ~ 10%, so that o Background — 10~ 4 — 10~ 2 fb, and it is at the level of 
the signal. However the kinematical configuration of the final state leptons is completely 
different. Figure [l(] (upper-right) shows the angular distribution of the gauge bosons which 
is peaked in the forward and backward directions so that the leptons produced by the W 
gauge boson decay are emitted preferentially along the collsion axis. In addition their 
tansverse momenta will be softer compared to that of the signal: Fig. |l0| (bottom-left 
panel) shows that the transverse momenta distribution of the gauge bosons is peaked at 
Pt = Mjv)/2 ~ 35 GeV for y/s = 300 GeV. Consequentely the lepton distributions 

will be peaked at ~ 17.5 GeV. The missing energy due to the undetected neutrinos 

(Fig. [l(], bottom-right panel) can be as large as ~ \/s — 2Mw- This distribution should 
be convoluted with that of the neutrinos produced in the gauge boson decay. Therefore it 
can be safely concluded that it will be possible to control this background because, with 
reasonable cuts on the transverse momenta of the leptons and on the missing energy, it will 
be drastically reduced, while - as mentioned above - these cuts will not affect significantly 
the signal. The same conclusion holds for the signal with heavy Major ana neutrinos: for 

> 700 GeV o~(2v2W) is 10 fb, the six particles final state has very large missing en- 
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Figure 10: Total cross section and distributions for e~e~ — ► W~W~vv. Upper-left figure: total 
cross section as a function of \fs. Upper-right: angular distribution for a W~ where 8 is the angle 
among the collision axis and the boson momentum. Bottom- left: distribution of the transverse 
momentum of W~ . Bottom-right: energy distribution of the two neutrinos. All distributions are 
calculated with V« = 300 GeV. 

ergy carried away by neutrinos and charged leptons have soft distributions in transverse 
momentum, thus cuts on missing energy and on the leptons pt reduce the background but 
not the signal. 

5. Conclusions 

The e~e~ option of the next generation of linear colliders offers the opportunity to test 
models of new physics through the discovery of lepton flavor violating signals, even if they 
arise only as a pure loop-level effect. We have shown, using the maximum experimentally 
allowed mixings, that masses of heavy Majorana neutrinos up to 2 — 3 TeV can be explored 
with the reaction e~e~ — ► £~£~, (£ = /tt,r), because the amplitude gets an enhancement at 
the treshold for two gauge bosons production and then shows a non-decoupling behaviour 
with the mass of the virtual heavy states. For the similar reaction e~e~ — > £~e~ , {£ = fi, r) 
induced by slepton mixing in supersymmetric models, in certain regions of the parameter 
space, the signal can reach the level of 10 -2 fb around the threshold for selectrons pair pro- 
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duction. The possibility of employing beams with high degree of longitudinal polarization 
is also essential to enhance the signal. On the other hand the standard model background 
is low and can be easily controlled. 
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